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Abstract 

The measurement of large scale anisotropies in cosmic ray arrival directions at energies above 10^'^ eV is 
performed through the detection of Extensive Air Showers produced by cosmic ray interactions in the atmo- 
sphere. The observed anisotropies are small, so accurate measurements require small statistical uncertainties, 
i.e. large datasets. These can be obtained by employing ground detector arrays with large extensions (from 
10^ to 10^ m2) and long operation time (up to 20 years). The control of such arrays is challenging and spu- 
rious variations in the counting rate due to instrumental effects {e.g. data taking interruptions or changes in 
the acceptance) and atmospheric effects {e.g. air temperature and pressure effects on EAS development) are 
usually present. These modulations must be corrected very precisely before performing standard anisotropy 
analyses, i.e. harmonic analysis of the counting rate versus local sidereal time. In this paper we discuss 
an alternative method to measure large scale anisotropies, the "East- West method", originally proposed by 
Nagashima in 1989. It is a differential method, as it is based on the analysis of the difference of the counting 
rates in the East and West directions. Besides explaining the principle, we present here its mathematical 
derivation, showing that the method is largely independent of experimental effects, that is, it does not require 
corrections for acceptance and/or for atmospheric effects. We explain the use of the method to derive the 
amplitude and phase of the anisotropy and we demonstrate its power under different conditions of detector 
operation. 
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1. Introduction 



The cosmic ray (CR) spectrum, in spite of its apparent regularity, exhibits in fact a few features, namely 
a first change of slope - the "knee" - near 10^^ eV, a second bending of the spectrum - the "second knee" - near 
10^^ eV, a hardening - "ankle" - around 10^^ eV, and a flux suppression around 5 x 10^^ eV. The measurement 
of the anisotropy in the arrival directions of cosmic rays is a complementary tool, with respect to the energy 
spectrum and mass composition, to investigate the origins of these features. From the observational point of 
view indeed, the study of the CR anisotropy, and especially its evolution over the energy spectrum, is closely 
connected to the problem of CR propagation and sources. 

Current experimental results show that the main features of the anisotropy are similar in the energy 
range (10^^ ^ 10^^ eV), both with respect to amplitude (lO""^ ^ IQ-^) and phase ((0 4) h LST) [l|,[3,ll,E 
EH- At higher ener gies , between 10^^ and 10^^ eV, the limited statistics do not allow any 
firm conclusion to be drawn \l3l . 14 , [la, , although the observation of a larger anisotropy amplitude 
with a difl'erent phase has been recently reported at ~ 4 • 10^^ eV (l7| . Around 10^® eV, when the gyro- 
radius of the particles becomes comparable to the galactic disk thickness, we would expect a large increase 
in the anisotropy {i.e. amplitudes at the percent level), in particular towards the galactic disk. From the 
experimental point of view, since the statistics are obviously more limited in this energy range, the situation 
is not as clear as at lower energies 18, [l^, [20I . 21 . 22|. Above 10^® eV, CRs (believed to be extra-galactic 



due to loss of confinement in the Galaxy) are not expected to show significant large scale modulations, with 
typical expectations for the amplitude arising from the Compton-Getting effect being below the one percent 
level. The change from a large to an almost null amplitude would mark the transition from galactic to 
extra-galactic CR origin. In few years the present experiments (e.g. the Pierre Auger Observatory (23| ) will 
have collected enough statistics to provide additional informations to probe this energy range. 

The experimental study of large scale CR anisotropics is thus fundamental for cosmic ray physics, though 
it is challenging. At energies > 10^'^ eV such measurements have to be indirect (due to too low primary 
fiuxes) and are usually performed through Extensive Air Shower (EAS) arrays. These arrays operate almost 
uniformly with respect to sidereal time thanks to the Earth's rotation: the zenith angle dependent shower 
detection and reconstruction is not a function of right ascension but it is a strong function of declination. 
Thus, the most commonly used technique (originally proposed in (23 |) is the analysis in right ascension only, 
through harmonic analysis (Rayleigh formalism) of the counting rate within the declination band defined 
by the detector field of view. Conventionally, one extracts the first and second harmonic: this is done by 
measuring the counting rate as a function of the local sidereal time (or right ascension) , and fitting the result 
to a sine wave. The Rayleigh formalism gives the amplitude of the different harmonics, the corresponding 
phase (right ascension of the maximum intensity) and the probability of detecting a signal due to fiuctuations 
of an isotropic distribution with an amplitude equal or larger than the observed one. 

The technique in itself is rather simple but the greatest difficulties are in the treatment of the data, i. e. 
of the counting rates themselves. Both for large scale anisotropics linked to diffusive motions or for the ones 
due to the Compton-Getting effect (i.e. due to the observer's motion with respect to a locally isotropic 
population of cosmic rays), the expected amplitudes are very small (10~^ — 10~^), with related statistical 
problems: long term observations and large collecting areas are needed. Instrumental effects must be kept as 
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small as possible, requiring detectors to operate uniformly (both in size and over time), and being as stable 
as possible. Moreover, EAS arrays are mostly located in remote sites (generally at mountain altitude), being 
thus subject to large atmospheric variations, both in temperature and in pressure. Meteorological induced 
modulations can affect the CR rate: indeed the EAS properties themselves depend on air density (through 
variations of the Moliere radius) and on pressure (due to the absorption of the electromagnetic component 
in the air) |25| . 

The measurement in practice is in fact complicated by the need of correcting the counting rate for 
instrumental and atmospheric effects, that must be done with high precision to prevent the introduction of 
artificial variations in the CR flux. 

The East-West method, being based on a differential technique, was designed to avoid introducing such 
corrections, preventing the possible associated systematics to affect the results. The original idea was pro- 
posed in [l|| and applied to the data of the Mt. Norikura array. It was later applied by other EAS arrays 
such as the Tibet experiment 0|, EAS-TOP 11| and the Pierre Auger Observatory 21, 22|. A modification 
of this technique was employed by the Milagro experiment [lO|. We revisit here this method, with the idea of 
explaining its mathematical background and of showing its power when applied to EAS data under different 
conditions of operation. In section [2] we explain the principle and show that the classical implementation 
of the method is valid within certain approximations. In section [3] we derive the mathematical basis of the 
method, avoiding some approximations used in the classical implementation and meanwhile demonstrating 
that the East- West technique is largely independent of any instrumental/atmospheric effect. The derivation 
of the amplitude and phase of the anisotropy from the harmonic analysis of the differences in East and West 
directions is illustrated in section 31 Here we show also how to extract from the derived amplitude and phase 
those corresponding to the equatorial component of the dipole. Before concluding in section [H we apply in 
section [5] the East- West method to different mock EAS data sets, characterized by different spurious effects. 



2. The principle of the East- West method and its classical implementation 

The East- West method relies on the fact that the difference between the observed counting rates of events 
recorded at each local sidereal time t (ranging from to 24 hs ; i.e. superimposing all detected events in a 
unique sidereal day) arriving from the Eastern and Western hemispheres, I'^^[t) and /^*(t) respectively, is 
proportional to the derivative of the true total counting rate ijot'^it), the coefficient of proportionality being 
approximately the mean hour angle (h) of the observed events. In this section, we aim at retrieving this 
classical implementation by outlining the different approximations which are needed to obtain this result. 
We will derive in the next section the relationship between I'^^{t) — I^^{t) and /j™^(t) in a more rigorous 
way. 

2.1. The principle 

The total counting rate of events observed in either the Eastern or the Western half of the field of view 
of an EAS array experiences different kind of variations during a sidereal day. Those may be caused either 
by experimental effects (changes of measurement conditions during the data taking, atmospheric effects on 
EAS, etc.) and/or by real variations in the primary CR fluxes from different parts of the sky. The East- West 
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method is aimed at reconstructing the equatorial component of a genuine large scale pattern by using only 
the difference of the counting rates of the Eastern and Western hemispheres. The effects of experimental 
origin, being independent of the incoming direction, are expected to be removed through the subtractior0. 
In the presence of a genuine dipolar distribution of CRs, as the Earth rotates Eastwards, the Eastern sky is 
closer to the dipole excess region for half a day each day; then, after the field of view has traversed the excess 
region, the Western sky becomes closer to the excess region and thus bears higher counting rates than the 
Eastern sky. The East- West differential counting rate is thus subject to oscillations whose amplitude and 
phase are expected to be related to those of the genuine large scale anisotropy. 

2.2. Classical implementation of the East- West method 

We denote the CR flux by $(a, (5) = + ^1(0., S), expressed in equatorial coordinates and where $0 
stands for the isotropic component, which is by construction large compared to the anisotropic one $1 (defined 
such that its average over the sky vanishes). The true counting rates (i.e. those that would be measured 
by a perfectly stable operating experiment and with negligible atmospheric effects on shower developments) 
from the Eastern and Western sectors, /|^"^(t) and /p|7"^(t), can be expressed in terms of the CR flux $(a, 5) 
as: 



where A is the effective area of the experiment and t stands hereafter, if not otherwise specified, for the 
local sidereal time and varies between and 27r. Ignoring spurious modulation effects that will be considered 
below, the experimental exposure uj{t,a,5) of a ground experiment as function of the local sidereal time t, 
the right ascension a and the declination 6 depends only on the combination (t — a, 6) in an even way with 
respect to the hour angle h = t — a. In the following, it will be useful to use the fact that once integrated 
over a in the field of view available at each time t, the resulting function uJ depends only on the declination: 



In real experiments, the observed counting rates I'^^{t) and I^^*(t) may be modulated by small instru- 
mental and atmospheric effects that influence measurement conditions during the data taking. Weather 
effects on EAS development actually modulate the estimation of the energy and thus the counting rate due 
to the steep energy spectrum, but they can formally be also treated as modulating in time the exposure of 



^This is a natural expectation provided the fact that the azimuthal detection efficiency of the corresponding experiment is 
symmetrical in East- West. Note however that corrections due to eventual asymmetries (such as those expected if the array is 
tilted) can be applied, as long as the asymmetries are known. 





(2) 




(3) 
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the experiment. For instance, if we adopt a simple modulation of the form u}{t — a, 5){l + r]{t)), with r]{t) the 
associated variation of the exposure, depending only on the experimental conditions at local sidereal time t 
and not on the direction {a, 6), the observed counting rates can thus be expressed as: 

ri'it) ~ I^i''%t)+A d6cos5 / dauj{t-a,5)'^or]{t), (4) 

r^max ft 

I^'{t) ~ d6cos6 / daLo{t-a,6)^ori{t), (5) 

where we neglected second order terms proportional to ^ir]{t). Since the exposure and the associated 
variations are identical in the East and West, it is straightforward to see that through the subtraction 
I'^^{t) — I^^{t), the terms proportional to r]{t) cancel, leading to: 



I'i'{t)-I^'{t) = A I d6 cos 6 



robs 



^max 



t+TT 

da uj{t — a, 5) 5) 



da u}{t — a, 5) $(a, 5) 



t-n 



(6) 



Assuming that the small angular scale variations of the flux "3>(a,(5) are small, we can approximate the 
integrations by considering the flux $(q,(5) as constant in each sector and evaluated at the right as- 
cension of the mean exposed direction at each time a = t i: (h), where (h) is the mean hour angle, 
(h) = Jl2"^ d6 cos 6 Jl_^ daio{t - a, 6){t - a). Then, 



/ d^cosJ 



^max 



^t + {h),S)-^t-{h),6) 



dauj{t — a, 5). 



(7) 



t-TT 



The additional simplification is to consider the mean hour angle small enough with respect to the scale of 
variation of so that one can use the linearised expressions <I>(t ± (h) , 6) ~ 6) it (h) d^{t, 5)/dt. Under 
these crude simplifications, the difference in the East and West counting rates can be approximated by: 



lt{t)-Iwit) - A I d5 COS 5 



jobs 



^max 



2{h) 



dt 



uJi6) 



^ {h) 



dt 



^max 



A I d6 cos6<^>{t,6)u}{6) 



(8) 



The expression between the brackets represents the true total counting rate, Itot^ii) = I^^^^{t) + I^r^^{t). 
Hence, we are led to the following relationship: 



(h) 



dt 



(9) 
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showing that the East- West counting rate difference, which is an observable quantity, is proportional to the 
derivative of the true intensity of cosmic rays, the coefficient of proportionality being approximately the 
mean hour angle, which is also a measurable quantity. 

To derive the classical result, it was necessary to perform some crude simplifications to get at Eqn.[8]from 
Eqn. O These approximations are expected to become less accurate as the zenithal range of the experiment 
is increased. However, we will see in the next section that the structure of the relationship still holds (with 
a different proportionality factor) in the general case. 



3. A more rigorous description of the East- West method 

We repeat here the previous calculations but without making the same simplifications as in the classical 
implementation. The aim is to check whether the relationship between the differential counting rate I'^^ — I^^ 
and the derivative of the total counting rate dl^^^^/dt still holds or not. 

It turns out that the calculation is easier when performed in local coordinates. The counting rates I'^^{t) 
and ly^'^ (t) at local sidereal time t for the two halves of the sky can be computed from the cosmic ray flux 
$ expressed in local coordinates [9, (p) as: 



/TT/Z /"fmax 
d0 / de smecos9e{9,t)^{9,(j),t), 
-n/2 Jo 

I^'{t) = A d(t) I d9 smecos9e{9,t)^{9,^,t), (10) 

Jtt/2 Jo 

where Acos9 is the effective area of the experiment at an angle of incidence 9 and e{9,t) is the detection 
efficiency function which includes the time-dependent spurious effects. Here, we adopt the convention that 
the azimuth angle is defined relative to the East direction, measured counterclockwise. To guarantee that 
the Eastern and Western sectors are equivalent in terms of counting rates, any dependence of e in azimuth (p 
needs to be symmetrical. For simplicity, we assume hereafter a uniform detection efficiency in azimuth; but 
similar conclusions still hold as long as the symmetry between the sectors is respected, which is a reasonable 
assumption in practice. It is also reasonable to assume that the relative amplitude r] of the temporal variations 
of the exposure is small, and that those variations decouple from the zenith angle dependent ones: 

e{9,t) = e^{9){l + v{t)). (11) 

To get an explicit expression for the cosmic ray flux <I> in local coordinates, we start from the param- 
eterisation $(a,5) in terms of the equatorial coordinates. The most basic approach to probe a large scale 
variation is to describe the flux by a combination of an isotropic component and a dipolar component: 

^{a,6) = %- (l + u{a,5)-D], D = D ■ u{ad,5d), (12) 



with D being the dipole vector deflned by its magnitude D and its orientation (ad, 5d) {u{a, 5) denotes 
the unit vector in the direction {a, 5)). As the conversion between equatorial and horizontal coordinates is 
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time-dependent, the flux from a specific viewing direction {9,(j)) expressed in local zenithal coordinates at 
any given location on the Earth turns out to be time-dependent. Adopting local coordinates with the z axis 
in the zenithal direction, the x axis towards the East and the y axis towards the North, the dipolar vector 
D{t) is written as 

/ D^{t) \ / sin 6d{t) cos Mt) \ 

D{t)= Dy{t) \ =D i smOdit) sin q)d{t) . (13) 
V D,{t) J \ cos 9d{t) J 

Transforming from local to equatorial coordinates, this can be expressed as: 

/ - cos 6d sin hdit) \ 

D{t) = D ■ cos I sin 5d — sinlcosSdCos hd{t) , (14) 
\ sin I sin 5d + cost cos 5d cos hd{t) J 

with hdit) = t — ad the hour angle of the dipole at the location of the observation point and at time t, and 
i the Earth latitude of the observation point. Inserting this expression in ^{6,(j),t) and carrying out the 
integrations over 6 and (p in Eqn. [10] leads to: 

/#^(i) = A^o-{^9n{t) + 2D^{t)gi2{t) + 7rD,{t)g2i{t)), (15) 
I^'it) = A<5o-{T^gn{t)-2D^{t)gi2{t) + TTD,{t)g2i{t)), (16) 

where the gij{t) coefficients are defined by: 



gij(t) = {l+rj{t))J dO ei{e) cos' esin> e (17) 
= [l+7l{t))hr (18) 

Once normalised by /n, each fij coefficient can be estimated directly from any data set through the empirical 
averages of cos*""*^ 6 ■ sin-'"^ 9: 

h. = /cos*-i e ■ sin-'-i e) . (19) 
/II 

As long as ^ixiax is smaller than — 70°, the zenithal distribution of CRs is only marginally affected by any 
large scale anisotropy patterns, in such a way that such empirical estimations, rigorous in case of isotropy, 
are still accurate enough even in the presence of genuine anisotropies. Then, by calculating as previously 
I'^^{t) — I^^{t) and dlf^^'^/dt, and by neglecting second order terms proportional to r]D, we find now that: 

ltit)-Iwit) ~ -4A<^oDh2Cos6dsinhd{t), (20) 



dt 

We then finally obtain: 



-2TTA^oDf2i cos £ cos 6d sin hd{t). (21) 



robsf,. robs... 2 (sing) dltont) 

Ie it) - % (t) - — (22) 
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Figure 1: Relative difFerences between the mean hour angle and the factor 2(sin^) /vr cos£ (cos as a function of the Earth 
latitude, for different conditions of maximal zenith angle available in the field of view. 



Consequently, the relationship between the East-West counting rate difference and the differential total 
counting rate is at first order similar to that discussed in the previous section, except that the proportionality 
factor is now given by /i^ = 2 {sin 9) / tt cos i {cos 9) , which can also be calculated from the measured zenith 
angles of the events. 

3.1. Comparison with the classical implementation 

Despite the crude approximations outlined in Section 12.21 it turns out that the factor /i^ is not very 
different from the mean hour angle as long as the Earth latitude of the site of the experiment is smaller than 
50° in absolute value. This is illustrated in Fig. [H where the relative difference {{h) — hi,)/hi, is shown as a 
function of the Earth latitude of the site, using e{9,t) = 1 and two different values of 9max- It can be seen 
that these relative differences become larger than 10% only for Earth latitudes i > 50°. The results would 
not be very different had we considered other examples of e{9,t). 



4. First Harmonic Analysis - Estimation of the dipole equatorial component 

J^.l. First Harmonic Analysis of Ie — Iw 

To probe a dipolar structure of the CR arrival direction distribution, Eqn. [22] is an ideal starting point 
to estimate the dipolar modulation of dll'^f^ /dt parameterised through the amplitude r and the phase ip: 

^l!^^ = rcos{t-ip). (23) 



8 



From a set of N arrival times from events coming from either the Eastern or the Western directions, r and 



(f can be estimated by applying to the arrival times of the events the standard first harmonic analysis [24 1 
slightly modified to account for the subtraction of the Western sector to the Eastern one. The Fourier 
coefficients a and b are thus defined by: 



2 ^ 



N 



2 



cos (ti + Q), (24) 



i=l 

N 



b = -Y^sm{U + Q), (25) 
1=1 

where Ci equals if the event is coming from the East or vr if coming from the WeslH. The amplitude and 
phase estimates {f,ip) of dlf^^'^/dt are then obtained through: 



vr cos i (cos 9) , ^ . f b\ 

-. T^y a + 0^ and ld = arctan - . 

2 (smy) \a J 



(26) 



By integrating dl/dt, the amplitude and phase estimates (f/, (pi) of the intensity I{t){^'^ itself are obtained: 

iV , vr 

r/ = — r and ipj =Lp + -. (27) 

ZTT Z 

Estimation of the dipole equatorial component 
In case of the standard Rayleigh analysis in right-ascension, the first harmonic amplitude r is related to 
the dipole amplitude D through [26|: 

(cos 5) Dj 



1 + (smO) D\\ 

where D\\ = D sinJ^ denotes the component of the dipole along the Earth rotation axis while D^^ = D cos 5^ 
is the component in the equatorial plane. The first harmonic amplitude thus depends on the declination of 
the dipole in such a way that it vanishes for 5^ = ib7r/2. This is obvious, as the modulation of the flux does 
not depend on the right ascension in such a case. On the other hand, the power of the method is the largest 
when the dipole is oriented in the equatorial plane. In this latter case the first harmonic amplitude becomes 
tra = D± {cos5) and the sensitivity of an experiment to the true value of D± depends then on (cos 5) , 
which is a function of the Earth latitude £ of the experiment and its detection efficiency in the zenithal range 
considered. 

Similarly, the first harmonic amplitude reconstructed by the East- West method is not directly the dipole 
amplitude. The additional step consists in relating dll^f^ /dt and the dipolar parameters through mea- 
surable quantities. Remembering that the observed number of events is equal to dt{I'^^ + I^r') = 



■^We are grateful to Paul Sommers for suggesting this simple way of accounting for the difference between the contributions 
from the East and West sectors. 
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27rA$o j2-Kd't^9iiif) + Dz{t)92i{t)) and that Dz{t) = DcosO^it) = D {sin i sin 6d + cos £ cos 5d cos hdit)), and 
neglecting all second order terms in riD± and ??^||, Eqn. [21] can be rewritten as: 

dP^rit) ^ iV D^(cosg)cos^sin/i^(0 
dt 2Trl + ^ J dtr]{t)+Di\ {cos e) sin i 

This expression describes how a genuine dipole with amplitude D and pointing to a declination 6d modulates 
in time the intensity of CRs observed by a ground experiment located at an Earth latitude i. The amplitude 
of the time variation is suppressed if the dipole is approximately aligned with the Earth rotation axis, while 
the modulation is maximal if is small compared to Dj^. In this latter case the differential counting rate 
can thus be expressed as a function of the equatorial component of the dipole in a straightforward way: 

Itit) - Iwit) ^ -^'^^^^D^ sinhdit), (30) 
zvr vr 

where we have neglected the second order terms proportional to r]D± and D^^D±. 

The first harmonic amplitude and phase estimates (D±,ad) are thus obtained through: 

= cos/{cos9) = nk9)'^'^^" «^ = ^ + f- (31) 

It is worth noting that, from the transformation of coordinates relation cos 5 sin /i = — sin ^ cos 93, the factor 
2 (sin^) /tt can be expressed as well in terms of 6 and h as: 

2 (sin 6) , ^ , , , . 
= (cos S sin h) , (32) 

IT 

where the r.h.s. average is performed by integrating over the eastern and western quadrants. 

By comparing Eqn. [31] to Eqn. [28] and using Eqn. [32l it can be seen that the use of the local sidereal time 
in the modulation search (instead of the right ascension as in the case of the standard Rayleigh analysis), 
combined to the East- West subtraction, leads to a loss of sensitivity by a factor (cos S) / (cos 6 sin h) (which 
is typically about two, depending on the experimental conditions) with respect to the performances of the 
standard Rayleigh analysis. However, this method has the benefit of avoiding the need to implement any 
corrections of the total counting rates for instrumental and atmospheric effects. Moreover, in some cases 
those corrections cannot be computed reliably, for instance when they are due to the energy dependence 
of the trigger efficiency, which can also depend on the unknown composition of the primary CRs. In these 
cases, only the East- West method can be implemented reliably. 



5. Simulations 

In this section, we check the behavior of the method through simulations reproducing realistic conditions 
of a ground experiment subject to artificial modulations at both the diurnal and the seasonal time scales. 
For this purpose, we consider three typical cases of interest for this kind of studies: 
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• isotropy; 

• genuine sidereal modulation; 

• genuine solar modulation. 

An additional underlying distribution has been considered to study the impact of a broken East- West sym- 
metry: 

• isotropy with a broken East- West symmetry. 

In order to test extreme situations we also add the spurious effects to all of these distributions. 

For definiteness we consider an experiment located at the same Earth latitude i = —35.25° as the 
Pierre Auger Observatory [2^ and sensitive up to a maximal zenith angle of 60°, with the following energy 
independent detection efficiency function: 

^) = T- n w T • 9{t), (33) 

l+exp{{e - 6ns)/(Tns) 

where 9ns = 50° and ans = 5°. Here, g{t) stands for the spurious modulation, which may differ from unity 
due to various reasons. For instance, the changes of atmospheric conditions affect the energy estimate of the 



showers [25|. When such effects are not accounted for, they induce a modulation of the rate of events above 



any given energy threshold. Hence, we choose the function g{t) to be of the generic form [2 



l + 7]y cos{2n{t -t^y)/Ty)+r]d 
9{t) = 



l+r/^cos(27r(t-tO)/Tj^) 



cos(27r(t - tO)/Td) 
, (34) 



where t is here expressed in terms of solar time. The cosine proportional to r]y describes the annual variation 
of the mean event rate, with phase and a period Ty of one year. The cosine proportional to rjd describes 
the annual average of the solar diurnal modulation, with phase t[] and a period of one day. Finally, r/y, 
stands for the variation of the diurnal amplitude along the year. This last term, combining the diurnal 
modulation with an annual one, is responsible for the production of sidebands at both the sidereal and the 



anti-sidereal frequencies, whose amplitudes are given by 0.5 x ry^ x ry^ [27|. To illustrate the power of the 
East- West method, we choose an extremely high value of 77^ = 90% to guarantee the existence of significant 
sidebands, while we set rjy = 20% and r/^ = 15%. 

5.1. Case 1: isotropy with spurious effects 

We first consider an isotropic distribution of CRs polluted by the spurious effects described above, and 
analyse 10^ mock samples generated with a total number of events N = 10^. The net results of the instru- 
mental and atmospheric effects introduced by the function g{t) are shown on top of Fig. [21 evidencing that 
the net counting rate = I'^'^ + undergoes the expected modulation of amplitude rjd = 15% at the 
solar frequency and of 0.5 x 77^ x r/^, = 6.75% at both the sidereal and the anti-sidereal frequencies. Applying 
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(e) 
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(f) 



Figure 2: Top: Distributions of the amplitudes of the first harmonic obtained with the standard Rayleigh analysis in the presence 
of experimental effects in solar time (a) and in both sidereal and anti-sidereal times (b). Middle and Bottom: Distributions of 
the amplitudes (c and d) and phases (e and /) of the first harmonic obtained with the East-West analysis, compared with the 
isotropic expectations: the subtraction of the spurious effects holds perfectly. 




Figure 3: Distributions of the amplitudes (a) and phases (fe) of the first harmonic in local sidereal time obtained with the 
East-West analysis in the presence of a genuine pattern of 5%. Both histograms are in perfect agreement with respect to the 
expectations. 



the East-West analysis at both the solar and sidereal frequencies, it can be seen that the reconstructed 
amplitudes Z)_|_ and phases are now distributed according to the expected distributions: the Rayleigh one 
for the amplitude with parameter a = 7r/2 (sin 6*) y^l/N, and the uniform one for the phase. Hence, in spite 
of the strong experimental effects, it turns out that the East- West subtraction allows the removal of possible 
biases in the estimate of both the amplitude and phase in the case of an underlying isotropic distribution of 
CRs. 



5.2. Case 2: 5% sidereal signal with spurious effects on the acceptance 

To test now the accuracy of the method in the presence of both a genuine signal at the sidereal frequency 
and spurious effects, a signal corresponding to a dipolar anisotropy of 5% amplitude is introduced in the 
simulated samples. For definiteness, we consider the dipole oriented towards the equatorial plane. The 
reconstructed amplitudes are now expected to follow a Rice distribution with parameters ji = 5% and 
a = 7r/2(sin6') yt/TV: 



M«x) = ^exp(^-^^J/„^-j^j, (35) 



while the reconstructed phases are expected to follow the distribution described by Linsley in the 2"^^ alter- 



native in 



1 / /^^ A , Mcos(arf-arf)/ ( ^1 cos {did - ad)W ( sin^ (a^ - a^) 
= ^ exp ^ - ^ J + 2V2^^ + ^^^ J J 



(36) 
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Figure 4: Same as Fig. |3] but generating a genuine pattern at the solar frequency. 



The results of the simulations are shown in Fig. [3l evidencing a perfect agreement with the expectations. 

5.3. Case 3: 0.1% solar signal with spurious effects 

We repeat here exactly the same exercise as above, but generating a genuine dipole with an amplitude 
0.1% at the solar frequency, together with the spurious effects. This kind of feature is expected due to the 
motion of the terrestrial observer through the frame in which the CR distribution is isotropic. It has been 
observed by several experiments at low energies, where sufficient statistics has been gathered. To probe such 
a low amplitude, the number of events has to be greatly increased with respect to the previous cases. Thus, 
we generated 10^ samples of = 10^ events. The results of the simulations are shown in Fig. HI showing 
once again perfect agreement with the expectations, even if both the genuine and the artificial modulations 
are present at the same time scale. 

5.4. Case 4-' isotropy with spurious effects and a broken East-West symmetry 

Finally, we study in this sub-section the influence of a broken East- West symmetry of the detection 
efficiency of a ground array experiment. There are several reasons in practice why the East- West symmetry 
may be slightly broken. We consider here an array being slightly tilted, the tilt being specified by a unit 
vector n described by its zenith angle 7 (the tilt) and its azimuth angle <j)Q (the direction). If the unit vector 
d describes a given direction in the sky with zenith 6 and azimuth (p the directional detection efficiency 
towards that direction is proportional to d ■ n: 

d - h = cos 7 cos + sin 7 sin 6 cos{(j) — c/^o) — cos 6* [1 + 7 tan 9 cos{<j) — cpo)] . (37) 
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Figure 5: Distribution of the amplitudes (a) and phases (6) of the first harmonic obtained with the East-West analysis in case of 
isotropy and of a broken East- West symmetry of the detection efficiency of the experiment. Thick histogram: sidereal frequency. 
Thin histogram: anti-sidereal frequency. Filled histogram: solar frequency. 



Introducing this azimuthal dependence into the detection efficiency results in a small difference between 
the Eastward and Westward counting rates, proportional to 'jcoscpQ at first order in such a way that for 
any North/South asymmetry the effect cancels exactly. Meanwhile, being independent of time, this shift 
does not impact itself in the estimate of the first harmonic. However, it is worth examining the effect 
of the combination of the tilted array together with the modulations induced by weather effects on EAS 
developments, because this combination leads to an East- West counting rate proportional to -yr]{t) cos <j)o, 
which may mimic a real East- West first harmonic modulation at the solar frequency. 

Generating 10^ isotropic samples of = 10^ events on a tilted ground array with a large value 7 = 5° 
and in the direction (j)Q = 0° (leading to the maximal effect), we show in Fig. [5] the results obtained at the 
sidereal frequency (thick histogram), the anti-sidereal frequency (thin histogram), and the solar frequency 
(filled histogram). It is clear that in such a case, only the solar frequency is affected by the East- West 
asymmetry introduced by the tilted array. On the contrary, the analyses performed at the sidereal frequency 
remains unbiased. 



6. Conclusions 

The differential East- West method for the measurement of large scale anisotropics has been revisited. 
Using the fact that the experimental instabilities simultaneously affect both the East and the West sectors, 
we have shown that with this method the equatorial component of the dipole can be recovered in an unbiased 
way, without applying any corrections to account for spurious effects, but with a reduced sensitivity with 
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respect to the standard Rayleigh analysis. Despite of this reduced sensitivity, this method has the advantage 
of avoiding the need to correct the total counting rate for instrumental and atmospheric effects. Finally, we 
have also shown that this method leads to unbiased results at the sidereal frequency even in the case of a 
broken East- West symmetry. 

This article is devoted to the memory of Gianni Navarra. Besides being a great scientist, he was deeply 
involved in this analysis and was the real inspirer of this work. 
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